Tidally-induced thermonuclear Supernovae 
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Abstract. We discuss the results of 3D simulations of tidal disruptions of white dwarfs by 
moderate-mass black holes as they may exist in the cores of globular clusters or dwarf galaxies. 
Our simulations follow self-consistently the hydrodynamic and nuclear evolution from the initial 
parabolic orbit over the disruption to the build-up of an accretion disk around the black hole. 
For strong enough encounters (pericentre distances smaller than about 1/3 of the tidal radius) 
the tidal compression is reversed by a shock and finally results in a thermonuclear explosion. 
These explosions are not restricted to progenitor masses close to the Chandrasekhar limit, we 
find exploding examples throughout the whole white dwarf mass range. There is, however, a 
restriction on the masses of the involved black holes: black holes more massive than 2 x lO'^ Mq 
swallow a typical 0.6 Mq white dwarf before their tidal forces can overwhelm the star's self- 
gravity. Therefore, this mechanism is characteristic for black holes of moderate masses. The 
material that remains bound to the black hole settles into an accretion disk and produces an X- 
ray flare close to the Eddington limit of LEdd — IC'^erg/s (Mbh/lOOOMa), typically lasting for 
a few months. The combination of a peculiar thermonuclear supernova together with an X-ray 
flare thus whistle-blows the existence of such moderate-mass black holes. The next generation 
of wide field space-based instruments should be able to detect such events. 



1. Introduction 

Two classes of black holes are established beyond reasonable doubt: black holes of a few solar 
masses ("stellar-mass") make themselves known by interacting with a companion star [1^ and 
" supermassive" black holes with masses of 10^ — 10^*^ Mq seem to lurk in the centres of most, if 
not all galaxies, see e.g. [2l|3]. A third class, so-called " intermediate- mass black holes" (IMBH) 
with masses of ~ 1000 Mq, represents a plausible, yet to date still missing link. The evidence for 
their existence is mounting, but to date it is not entirely conclusive. Based on kinematical studies 
Gerssen et al. [H |5] have suggested the presence of a 3.9 x 10^ Mq black hole in the core of the 
globular cluster M15. More recently, there were further reports on IMBHs in the globular clusters 
NGC 2808 [6j and omegaCen [7j. Additional evidence comes from ultraluminous, compact X-ray 
sources in young star clusters [H [9] and from n-body simulations [10] that indicate that runaway 
collisions in dense young star clusters can produce rapidly growing massive stars that should 
ultimately collapse to form IMBHs. None of the above arguments is strong enough to close the 
case, but their different nature gives this hypothesis credibility. 

At the very large number densities in the environments that may habor such black holes 
stars are disrupted at a rate of ~ 10~''yr~^ ^hCs ^*:6 crio (.^per/Rt), where Mbh,3 is the 
black hole mass in units of 1000 Mq, n^K^g is the central stellar density in units of 10^ pc~^, 



(Tio the velocity dispersion in 10 kms and -Rper is the distance of closest approach. The 
quantity i?t = (Mbh/M,,,)^''^ i?,,, is the tidal radius, the approximate separation where the 
hole's tidal field wins the battle against the star's self-gravity. It is proportional to the 

1/3 

stellar radius i?*, but only grows oc M^^ , i.e. substantially slower than the gravitational 
radius, Rs = 2GMbh/c^, of the black hole. As a consequence, too massive black holes, 

cr 3/2 —1/2 

Mbh > Mbh,iim — 2.5 X 10 Mq ^„d9^^wd 6 ^^^^ swallow stars before disrupting them. Here, 
Rwd,9 is the white dwarf radius in 10^ cm and -Mwd,o.6 its mass in units of 0.6 Mq. Thus, white 
dwarfs with their small radii and their (at least in principle) available nuclear fuel can only be 
tidally squeezed and disrupted by black holes of moderate masses. Therefore they are precious 
tools to probe the existence of this type of black hole. 

Luminet and Pichon were the first to discuss the tidal disruption of white dwarfs by black 
holes [11] in the context of the affine star model [12| [13]. Wilson and Mathews [l^ believe 
that general relativistic effects may lead to a tidal ignition even for supermassive black holes 
and at separations of many gravitational radii. A hydrodynamic simulation with a Newtonian 
Adaptive Lagrangian Eulerian (ALE) code was performed recently [15] where the gravitational 
constant was adjusted to mimic general relativistic effects. We have recently performed a large 
set of 3D hydrodynamic simulations |16l [T7| [T8] that incorporate the nuclear energy generation 
self-consistently with the hydrodynamics. 

2. Ingredients for the simulations 

Here we briefly summarize the ingredients of our simulations, both in terms of physics and 
numerical methods. We described the numerical aspects in [T7] and in [18] we concentrated on 
the physics and the astrophysical implications of such tidal disruptions. 

We use the smoothed particle hydrodynamics method (SPH) to follow the hydrodynamic 
evolution of the white dwarf matter. Being completely Lagrangian and conserving mass, energy, 
momentum and angular momentum by construction the method is ideal for the study of such 
a violent stellar disruption. General reviews on the method can be found in \19\ [20] [2T] , our 
particular implementation is documented in [T7]. We calculate the forces from self-gravity via 
a binary tree |22] and those from the central black hole via a relativistic pseudo potential [23], 
details of how the singularity is treated numerically are laid out in [2?! .We close the system 
of hydrodynamics equations via the HELMHOLTZ equation of state (EOS) ^25j. It accepts 
an externally calculated chemical composition, facilitating the coupling to nuclear reaction 
networks. The electron-/positron equation of state makes no assumptions about the degree of 
degeneracy or relativity, the exact expressions are integrated numerically and are subsequently 
tabulated. The interpolation in this table enforces thermodynamic consistency by construction 
[25]. The nuclei in the gas are treated as a Maxwell-Boltzmann gas, the photons as blackbody 
radiation. Most importantly in this context, we couple a minimal nuclear reaction network [26] 
to the hydrodyamics. Throughout this simulation set we assume a uniform nuclear composition 
accross the white dwarfs. Stars with masses < 0.6 Mq are instantiated as pure helium, more 
massive stars are modeled with a mass fraction of 50% carbon and 50% oxygen. All simulations 
are started from parabolic orbits, the strength of an encounter is parametrized by the penetration 
factor (3 = i?t/-Rper- For more details we refer the interested reader to |tl7( il8j . 

3. Results 

3.1. Explosion mechanism 

The tidal gravitational field is, to lowest order, determined by the second-order tidal tensor 
and the third-order deviation tensor. The latter determines the deviation of an extended body 
from a point particle trajectory of the same mass. The tidal tensor determines how the fluid 
will be shaped. It possesses one positive and two negative eigenvalues, see e.g. [2T, and, as 



a consequence, the tidal field stretches the star in one eigen- direction and compresses it along 
the other two. Most severely, the star is compressed perpendicular to the orbital plane, this 
process halts when a shock forms \2S\ ?I7\ [TH] that reverts the compression into an expansion. 
For deep enough penetrations matter reaches nuclear statistical equilibrium, see Fig. [T] upper 
right panel. The star is squeezed through a point of maximum compression on a time scale 

of Tcomp ~ R^d/Vp ~ 0.2s (M^d,0.6)"'/'(^wd,9)'/'(Mbh,3)~'/', where Vp ~ {RjRr)'/^c ~ 

5 X lO^cm s~^ -^wd^oe ^wdd^hCs peri-centre passage velocity. This velocity exceeds 

thermonuclear flame speeds by orders of magnitude, therefore, flame propagation effects can be 
safely neglected. The compression time scale needs to be compared to the dynamical time scale 
of the star, Tdyn = (Gp)^^/^ ~ 7.2 s M^^qq Rwdg^ with p being the average stellar density 
and to the nuclear reaction time scale, Tnuc- Only for Tquc *C Tcomp, Tdyn can a substantial 
nuclear energy release be expected. Our large simulation set |18] shows that white dwarfs of all 
masses can be thermonuclear ly exploded provided that the penetration factor exceeds /3 ~ 3. 
For definiteness, we illustrate here the results at the example of a typical 0.6 M0 carbon-oxygen 
(CO) white dwarf that passes a 500 M0 black hole with a penetration factor of /3 = 5, see Fig. 

m 

3.2. Nucleosynthesis 

Nucleosynthesis is triggered predominantly in the point of maximum compression, see Fig. [T] 
upper right panel, the concomitant nuclear energy release inflates a hot bubble in the debris 
centre, see Fig. [T] lower panels. The high-density centre of the star produces an iron-group core 
of 0.18 if mainly composed of ^^Ni this amount is comparable, but at the lower end of 

what is deduced for standard type la supernovae [29]. The iron core is surrounded by a 0.21 Mq 
shell of silicon-group elements which in turn is covered by a sheath of unburned carbon-oxygen 
material, see Fig. [2] In reality, such a disruption would result in even stronger carbon-enhanced 
outer layers than shown in Fig. [2j For simplicity, we have instantiated our initial carbon oxygen 
white dwarf models (Mwd ^ 0.6 Mq) as homogeneously mixed stars with a 50% mass fraction 
of each nucleus. While such internal chemical profiles are likely accurately realized in nature 
in very massive white dwarfs (~ 1 M©) [30], for lower masses the gravothermal adjustment 
of the interior during the cooling phase produces oxygen-enhanced stellar cores surrounded by 
very carbon-rich mantles {Xc ~ 0.8). The exact radial distribution depends on the exact value 
of ^^C(a,7)^^0 rate and the details of how convection proceeds |30| [3T| I32j . but this general 
stratification tendency is well-established. Thus, the disruption of a standard 0.6 M0 white 
dwarf should produce a highly carbon-enriched remnant atmosphere. 

We found exploding examples throughout the whole white dwarf mass range. An extreme case is 
a low-mass white dwarf that can (due to its large stellar and correspondingly large tidal radius) 
be sent very deeply into the tidal radius of a black hole. A 0.2 Mq, pure helium white dwarf 
disrupted by a 1000 Mq black hole (/3 = 12) produced an explosion with about 0.03 Mq of 
iron-group elements. At the other extreme, a 1.2 Mq CO white dwarf exploded with 0.66 Mq 
of iron group elements after passing a 500 Mq black hole with (3 = 2.6. For more details and 
examples, the interested reader should consult [18j. A detailed comparison of the lightcurves of 
tidally-induced supernovae and "normal" type la supernovae is currently being prepared [33] . 

3.3. Accompanying signatures 

If the release of nuclear binding energy is disregarded, about half the stellar mass is ejected, 
the rest is still bound to the black hole and, in principle, available to be accreted. Nuclear 
energy release during pericentre passage can substantially increase the amount of unbound mass 



^ Note that our network uses element groups, individual isotopes would need to be recovered in a post-processing 
step with larger network. 




Figure 1. Density (left column) and temperature (right column) evolution during the disruption 
of 0.6 M0 CO white dwarf by a 500 M0 black hole. The distance of closest approach was 1/5 
of the tidal radius. 




Figure 2. Mass fractions resulting from the disruption of 0.6 Mq CO white dwarf by a 500 Mq 
black hole. The distance of closest approach was 1/5 of the tidal radius. 



to about ~ 65% [16J. Before the energy contained in the bound matter can be released, an 
accretion disk has to form. On returning to the black hole, the orbits of the bound matter 
become focussed once more towards peri-center. The spread of specific energies across the 
accretion stream is very large [34j so that matter that passes on the inside of the stream, i.e. 
closer to the black hole, is considerably stronger bound and subsequently only reaches moderate 
apocentre distances while matter on the outside track is relaunched into a substantially wider 
orbit. Thus, after having passed the hole for the second time, the matter stream spreads in a 
fan-like manner |34| [T8] . It subsequently collides with the still infalling material and this self- 
interaction produces an angular momentum redistribution shock that circularizes the accretion 
flow, see e.g. Fig. 4 in [T7]. In the early stages, mass is fed towards the black hole at a rate that 
carries the imprint of the internal structure of the star |35l |36], at later stages the rate falls off 
oc t~^/^ [Ml EZ]. Once a disk has formed, it evolves under the influence of viscosity, radiatively 
driven winds and the fallback rate by which it is fed. We expect a luminosity that is comparable 
to the Eddington value, LEdd — lO^^erg/s (Mbh/lOOOMo), typically lasting for a few months. 

4. Summary and discussion 

We have explored in detail the fate of white dwarfs that approach a moderately-massive black 
hole close enough to become tidally disrupted. Above a limiting mass of ~ 10^ Mq, the detailed 
value depending on the white dwarf, the black hole swallows the star as a whole before disrupt- 
ing it. Thus white dwarfs represent a precious tool to probe the existence of moderate-mass 
black holes. White dwarfs that penetrate the tidal radius deeply enough are heated by the 
tidal compression and ensuing shock to the temperatures of nuclear statistical equilibrium. In 
favorable cases more than the gravitational self-binding energy of the star can be released via 
nuclear reactions thus triggering a thermonuclear explosion of the white dwarf. The amount of 
iron-goup nuclei produced is rather sensitive to the densities at which nuclear burning occurs, 
and therefore to the white dwarf mass. The amount produced in the exploding cases of our 
simulation set |18] ranges from 0.03 Mq for a 0.2 Mq pure helium white dwarf to 0.66 Mq for a 
1.2 Mq carbon-oxygen white dwarf, i.e. if composed of mainly ^^Ni the favorable cases produce 
amounts comparable to standard type la supernovae [29j. 

These tidally-induced supernovae are, however, in several respects different from what is con- 
sidered a "normal" type la supernova. First, they are not restricted to progenitor masses close 
to the Chandrasekhar limit. We find exploding examples throughout the full explored mass 
range from 0.2 to 1.2 Mq. Typical examples of these explosions are close to the peak of the 
white dwarf mass distribution, i.e. near 0.6 Mq, or slightly above this value, since in the dense 
stellar environments around such black holes close encounters between stars may lead to mass 
segregation [38j. Thus the average progenitor will be less massive than for a normal type la. 
Second, and closely related to the first point, the progenitor composition is not restricted to 
carbon-oxygen. Also pure helium white dwarfs can explode via this mechanism, but since white 
dwarfs near and slightly beyond the peak of the mass distribution are thought to consist of 
carbon and oxygen this is the most common progenitor composition. Due to the initial chemical 
layering of the white dwarfs, the outer shells of the most common explosion remnants consist of 
unburned, highly carbon-enriched [Xq ~ 0.8) material. Third, the remnant geometry is most 
likely much less spherical than a standard type la, therefore the optical lightcurve should be 
rather unique as a result of the radiating material being highly squeezed into the orbital plane. 
The large velocities of the ejected debris (> 10^ km/s) should produce large Doppler shifts. And 
last but not least, this peculiar type of thermonuclear supernova is accompanied by an X-ray 
flare close to the Eddington luminosity that lasts for a few months. 

The estimated event rates for this type of transient are ~ 10^'^ of the type la supernova rate. 
Though being substantially less frequent, they occur often enough to warrant a search for this 
new class of optical transient. Upcoming supernova searches hope to discover several thousand 



to several hundred thousand type la-Hke events |39| HO] per year. Chances are promising to 
find examples of tidally-induced thermonuclear supernovae among them. Maybe, the recently 
detected carbon-rich transient SCP 06F6 accompanied by an X-ray signal jJT] is already the 
first example for this class of object. 
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